Recent progress 
on nuclear parton distribution functions 



M. Hirai*, S. Kumano* and K. Saito* 

* Department of Physics, Faculty of Science and Technology, Tokyo University of Science 
2641, Yamazaki, Noda, Chiba, 278-8510, Japan 
' KEK Theory Center, Institute of Particle and Nuclear Studies, KEK 
and Department of Particle and Nuclear Studies, Graduate University for Advanced Studies 
1-1, Ooho, Tsukuba, Ibaraki, 305-0801, Japan 

Abstract. We report current status of global analyses on nuclear parton distribution functions 
(NPDFs). The optimum NPDFs are determined by analyzing high-energy nuclear reaction data. Due 
to limited experimental measurements, antiquark modifications have large uncertainties at x > 0.2 
and gluon modifications cannot be determined. A nuclear modification difference between u and 
d quark distributions could be an origin of the long-standing NuTeV sin 2 9 W anomaly. There is 
also an issue of nuclear modification differences between the structure functions of charged-lepton 
and neutrino reactions. Next, nuclear clustering effects are discussed in structure functions as a 
possible explanation for an anomalous result in the 9 Be nucleus at the Thomas Jefferson National 
Accelerator Facility (JLab). Last, tensor-polarized quark and antiquark distribution functions are 
extracted from HERMES data on the polarized structure function b\ of the deuteron, and they could 
be used for testing theoretical models and for proposing future experiments, for example, the one at 
JLab. Such measurements could open a new field of spin physics in spin-one hadrons. 

Keywords: Quark, gluon, parton, distribution, QCD, nuclear effect 
PACS: 13.60.Hb, 13.60.-r, 24.85.+p, 25.30.-c, 13.88.+e 



INTRODUCTION 

Parton distribution functions (PDFs) of the nucleon have been investigated for a long 
time, and they are now established except for extremely small-x and large- x regions. A 
precise determination of the PDFs is very important for new discoveries at LHC (Large 
Hadron Collider). For heavy-ion collisions at LHC, nuclear parton distribution functions 
(NPDFs) should be also determined precisely because 10—20% nuclear modifications 
exist in medium-size nuclei and they are larger in heavy nuclei. In addition to heavy-ion 
physics at LHC, an accurate determination of the NPDFs is valuable for other topics such 
as establishment of non-perturbative aspect of QCD by comparisons with theoretical 
models of nuclear medium effects, neutrino-oscillation studies by nuclear corrections in 
16 0, high-energy cosmic-ray interactions with air for studying GZK (Greisen-Zatsepin- 
Kuzmin) cutoff, and determination of fundamental constants like sin 2 Oyy- In contrast to 
the nucleonic PDFs, the NPDFs have not been determined well, which adds uncertainties 
in investigating the above-mentioned topics. In this report, we explain the current status 
of the NPDF determination and our recent studies on nuclear structure functions. 

Nuclear modifications of the structure function F2 have been investigated in deep in- 
elastic lepton scattering from nuclear targets [1]. Now, we have measurements from 
relatively small x to large x with a variety of nuclei from the deuteron to heavy ones. Us- 



ing these data together with Drell-Yan measurements of nuclear targets, we determined 
the NPDFs. RHIC (Relativistic Heavy Ion Collider) and neutrino reaction data could 
be used in addition. First, we introduce the current status of the NPDF determination 
[2, 3, 5, 4, 6, 7] and a possible relation to the NuTeV sin 2 Ow anomaly [8, 9, 10]. We 
also discuss an issue of different nuclear modifications between charged-lepton and neu- 
trino reactions. Second, the structure function Fi of 9 Be is discussed for explaining an 
anomalous JLab data on the nuclear modification slope d(F 2 Be /F 2 D )/dx [11]. Since 9 Be 
is a nucleus which is usually described by clustering structure of two a nuclei with neu- 
tron clouds, the data could indicate a clustering feature in the deep inelastic scattering 
(DIS) [12]. Third, new spin structure of the deuteron is investigated by analyzing HER- 
MES data on the tensor-structure function b\ [13, 14]. The tensor structure of nuclei and 
hadrons in the high-energy region is an unexplored topic, although it has been investi- 
gated extensively at low energies in connection with nuclear tensor force due to meson 
exchanges. It is interesting to understand the tensor structure in terms of quark and gluon 
degrees of freedom. Since b\ vanishes if constituents are in the 5-wave, it could also be 
a good probe of quark-gluon dynamics including orbital angular momenta. 



Nuclear modifications of the structure function F2 are typically 10—20% for medium 
size nuclei. Therefore, it is convenient to express the NPDFs by the corresponding 
nucleonic PDFs multiplied by nuclear modification factors w,- (z = u v , d v , q, g): 



where p and n indicate the proton and neutron, A, Z, and N are mass, atomic, and 
neutron numbers, respectively, and Qq is the initial Q 2 scale. The isospin symmetry 
(d n = u p , u n = d p ) is assumed in relating the distributions of the neutron to the ones of 
the proton. It should be noted that the above parametrization cannot describe the NPDFs 
at x > 1 because the distributions vanish in the nucleon, whereas the kinematical range 
is < x < A in a nucleus. 

There are a few groups which have been investigating the optimum NPDFs by global 
analyses of world data on high-energy nuclear reactions [2, 3, 5, 4, 6, 7]. The x- 
dependent functional forms are different depending on the analysis groups. For example, 
the HKNS parametrization uses [2] 



at Qq = 1 GeV 2 . Here, a, /3, a ; , bu Ci, and d[ are parameters to be determined by a 
analysis. Other groups' functional forms are listed in Ref. [15]. In the parametrization 
of de Florian and Sassot [3], the NPDFs are expressed by convolution of the nucleonic 
PDFs with a modification function w(x) and a small initial Q 2 scale (0.4 GeV 2 ) is used 
for utilizing the GRV-type PDFs. In the analysis of Eskola, Paukkunen, and Salgado 
(EPS) [4], the x region is divided into three: shadowing, anti- shadowing, and EMC & 



NUCLEAR PARTON DISTRIBUTION FUNCTIONS 



ff(x,Q 2 ) = Wi (x,A,Z)j [Zff(x,Q 2 )+Nf?(x,Q 2 )} , 



(1) 




Fermi-motion regions, and different parameters are assigned to the three functions in 
each NPDF. Schienbein et a/.'s analysis used a functional form of polynomials of x 
combined with exponential functions [5]. Used data sets are similar; however, the EPS 
[4] included some of the RHIC data and Schienbein et al. [5] also used NuTeV neutrino 
data. There are other works by Kulagin and Petti (KP) [6] and Frankfurt, Guzey, and 
Strikman (FGS) [7]. The KP method is to determine off-shell structure functions of the 
nucleon by calculating other nuclear modifications in conventional models with binding, 
Fermi motion, and multiple scattering. The FGS method is mainly intended to describe 
the small-* shadowing part by the multiple scattering. 



fi p s b 




FIGURE 1. Nuclear modifications of the PDFs are shown for the lead 
nucleus at g 2 =1.69 and 100 GeV 2 [4]. The functions R™, R P s h , and 

R^ b indicate the valence-quark, sea-quark, and gluon modifications, 
respectively. Three parametrization results (EPS09, HKN07, and nDS) 
are shown with uncertainty bands. 



Nuclear modifications of the PDFs are determined by % 2 analyses of world data on 

the structure-function F2 ratios F^ jF^ , Drell-Yan cross-section ratios ojjy/ G DY > RHIC 
dA data, and neutrino-DIS structure functions. Typical results are shown in Fig. 1 for 
nuclear modifications of the lead nucleus at Q -1.69 and 100 GeV 2 . Three different 
parametrizations are shown together with uncertainty bands. Although the three groups 
used different parametrizations for the x dependence and slightly different data sets, the 
determined distributions roughly agree with each other within the uncertainty bands. 
There are large uncertainties in the antiquark distributions at medium and large x and 
in the gluon distribution at whole x. They should be determined by future experimental 
measurements. 

The studies of nuclear modifications could affect the NuTeV anomaly on the weak- 
mixing angle. The NuTeV anomaly has been a long-standing issue that the mixing angle 
measured by the NuTeV neutrino DIS, sin 2 W = 0.2277 ± 0.0013 (stat) ± 0.0009 (syst), 
is different from other measurements, sin 2 &w = 0.2227 ±0.0004 [8]. In extracting the 
mixing angle from neutrino DIS, there is a Paschos-Wolfenstein (PW) relation for the 
isoscalar nucleon. There are various corrections to use it for the iron target in the NuTeV 



experiment: 



vA _ G vA j 

f A _ f A =■=- sin 2 w + (Z^N) + (w Uv ^ w dv ) + (u p ^ d n ) + (s± s) + (c ^ c) , (3) 
°cc °cc 

where and are charged-current (CC) and neutral-current (NC) cross sections, 
respectively, in neutrino-nucleus DIS. The parentheses (• • •) indicate various corrections 
to the PW relation. The obvious non-isoscalar correction (Z 7^ N) is taken into account 
in the NuTeV analysis, so that it should not be the origin of the anomaly. 

First, since the NuTeV target is the iron nucleus, the anomaly could come from 
nuclear modification difference between u v and d v . In particular, baryon-number and 
charge conservations indicate that the modifications should be different between u v 
and d v [9]. However, the current data are not enough to probe such a small effect 
by a global NPDF analysis [9]. According to Cloet, Bentz, and Thomas, the nuclear 
modification is larger for u quark than the one for d due to vector mean field which is 
adjusted to explain the symmetry energy of a nucleus [10]. Future pion-induced Drell- 
Yan measurements such as by the COMPASS collaboration could be used for finding 
the difference [16]. Second, QED effects on the PDFs, namely u p 7^ d n and d p 7^ u n 
could be the origin of the anomaly. The symmetry relations between the proton and 
neutron PDFs (u p = d n , d p = u n ) are usually assumed in all the global analyses by 
considering that the QED effects are less than a few percent (~ Oc/oc s ). A possible QCD 
effect was investigated in Ref. [17] by somewhat assuming a photon distribution at the 
initial scale Q 2 =l GeV 2 . The result indicated that the anomaly could come from such an 
effect although a precise determination of the QED effect is not possible at this stage. 
Third, the strange and charm asymmetries (s 7^ s, c 7^ c) could be origins. However, such 
effects cannot be determined from experimental measurements at this stage in the sense 
that the distribution s(x) — s(x) still has a large uncertainty band [18]. 

Another issue of the current NPDFs 
is the modification differences between 
charged-lepton and neutrino reactions [5]. 
The nuclear modifications of F2 are well 
established in charged-lepton reactions, 
and the optimum NPDFs are extracted 
mainly from their data. However, the 
modifications are very different in neu- 
trino DIS as shown by the solid curve in 
Fig. 2 [5], whereas the other curves are 
the modifications suggested by charged- 
lepton DIS. The medium-* (= 0.5 ~ 0.8) 
modifications are different and the anti- 
shadowing phenomenon does not appear 
at x = 0. 1 ~ 0.2 in the neutrino DIS. Since 
neutrino data are partially used for extracting current "nucleonic" PDFs, this issue should 
be solved for accurate determinations of not only the NPDFs but also the PDFs of the 
nucleon. 
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FIGURE 2. The solid curve (fit A2) indicates 
the nuclear modification of i<2(iron) in neu- 
trino DIS [5]. For comparison, two parametriza- 
tion results (KP and HKN07) are shown to- 
gether with the modification in charged-lepton 
DIS (SLAN/NMC). 



CLUSTERING EFFECTS IN NUCLEAR STRUCTURE 

FUNCTIONS 



It is known in nuclear- structure studies that some light nuclei exhibit phenomena of 
nuclear clustering. For example, the 8 Be nucleus is theoretically described by two a- 
like clusters according to a Monte Carlo calculation [19] rather than the usual monotonic 
density distribution estimated by a shell model. This fact led us to wonder whether such 
clustering structure could appear in high-energy reactions, for example, in the structure 
functions F 2 - Although simple multiquark clusters were investigated at the early stage 
of EMC-effect studies, clustering aspects of nucleons were not investigated. Recently, 
such a signature could be found by a JLab experiment [11] according to the theoretical 
work of Ref. [12]. We explain this theoretical approach on the clustering effects. 

First, we estimate effects of nuclear clustering by using a conventional convolution 
approach. A nuclear structure function F 2 is given by the nucleonic one convoluted 
with a nucleon momentum distribution in the nucleus f(y): [1, 20] 

Fi(x,Q 2 ) = jy y f(y)F 2 N (x/y,Q 2 ), f(y)=jjd'p N y8^y-^-^\HPN)\ 2 , (4) 

where F 2 is the structure function for the nucleon, and y is the momentum fraction y = 
Ma Pn ■ q/ {M N pA-q) ^Ap~^/ [p\) with a light-cone momentum p + . If the momentum 
density \(j)(p N )\ 2 is obtained, F 2 can be calculated by this equation. 

Two theoretical models are used for calculating the nuclear densities. In order to 
investigate a possible clustering effect, an antisymmetrized molecular dynamics (AMD) 
is used and a simple shell model is also employed to compare with AMD calculations. 
The AMD is a variational approach for describing nuclei. Its advantage is that there is no 
a priori assumption on nuclear structure whether it is a shell or cluster-like configuration. 
A nuclear wave function is given by the Slater determinant of single-particle wave 
packets: \®(r h r 2 ,- ■ -,r A )) = (l/v^)det[pi(ri), <^(r 2 ), • • •, 9a(?a)] where <p ; (r ; ) is the 
single-particle wave function given by <Pi(rj) = (j>i(rj) Zi T i w i m spin and isospin states Xi 
and ti. The function 0j(r/) is the space part expressed as 0/(r/) = (2v/7r) 3 / 4 exp[— v(r ; - — 
Zi/y/V) ], where v and Z,- are variational parameters. Simple effective AW interactions 
are used for describing gross properties of nuclei. The variational parameters are then 




FIGURE 3. Two-dimensional (x, y) coordinate- FIGURE 4. Two-dimensional (x, y) coordinate- 
space density of 4 He by the AMD [12]. The den- space density of 9 Be by the AMD [12]. The nota- 
sity is shown by taking an integral over the coor- tions are the same as the one in Fig. 3. The shell - 
dinate z: J dzp (x 7 y,z). The shell-model density is model density is a monotonic distribution, which 
almost the same as the above distribution in 4 He. is different from the above distribution. 



determined by minimizing the system energy with a frictional-cooling method. Obtained 
space densities are shown in Figs. 3 and 4 for 4 He and 9 Be. Figure 3 indicates that 
the 4 He density is a monotonic distribution as typically suggested by shell models. 
However, it is obvious from Fig. 4 that the 9 Be distribution is totally different. There 
are two density peaks which roughly correspond to two a clusters with surrounding 
neutron clouds. It is interesting to investigate that such a nuclear clustering could be 
found in high-energy experiments, for example, in nuclear structure functions because 
the anomalous EMC effect was indicated by the JLab experiment for 9 Be [1 1]. 

In order to show the clustering effects, a simple shell-model is used for calculating 
nuclear wave functions in comparison with the AMD results. We use a harmonic- 
oscillator model with the potential M^co 2 r 2 /2. Its wave function is given by i//„^ m = 

R n e(r)Y hn (0, 0), Rni(r) = A^rV^ r L n _{ (kV), where Y bn (0 : (j)) is the spherical 

harmonics, ( x ) * s tne Laguerre polynomial, K is defined by K = y/M^co, and Nni is 

the normalization constant. The shell-model density of 9 Be is a monotonic distribution, 
which is totally different from the AMD one in Fig. 4. However, if an angular average is 
taken for the density, the difference between the AMD and shell model is not very large 
[12]. Taking the Fourier transformation to the momentum space, we show both densities 
in Fig. 5 for 4 He and 9 Be. Although the AMD and shell-model densities are almost 
the same in the 4 He nucleus, they are different in 9 Be. The AMD distribution is shifted 
toward the high-momentum region, which is caused by the nuclear clustering. Since 
nucleons are confined mainly in two small space regions (clusters), the high-momentum 
components are created. It is interesting to find such a difference due to the nuclear 
clustering. We then estimate its effect on nuclear structure functions by using Eq.(4). 

The EMC ratio (F^ /F^) is shown for the 9 Be nucleus in Fig. 6 together with the 
JLab data. Since effects of short-range correlations and internal nucleon modifications 
are not included in our estimate, the obtained curves do not fully agree with the data. 
Here, it is our intention to show the clustering effect in the basic convolution formalism. 
Considering that the differences between curves are of the order of the experimental 
errors and there are some theoretical ambiguities, for example, from the correlations 




FIGURE 5. Nucleon-momentum distributions in FIGURE 6. Nuclear modifications of for 9 Be 

4 He and 9 Be by the shell and AMD models [12]. [12]. The dashed and solid curves indicate the 

They are obtained by taking Fourier transforma- shell-model and AMD results at Q 2 =5 GeV 2 . 

tions of the coordinate-space densities. JLab measurements are shown for comparison. 
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FIGURE 7. Nuclear modification slopes mea- 
sured at JLab are shown by taking the scaled nu- 
clear density as the abscissa [11]. The curve in- 
dicates a fit to the data except for 9 Be. 
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FIGURE 8. Nuclear modification slopes are 
shown by taking the maximum local densities in 
the shell and AMD models as the abscissa [12]. 
The curve indicates a fit to the data points of the 
shell model except for 9 Be. 



and nucleon modifications if they are included in our estimate, it may not seem easy to 
find a clear clustering signature. 

The Jlab experimental group obtained the nuclear modification slope dR/dx for 
R = F£ /F!? by approximating their data by straight lines at 0.35 < x < 0.7 [11]. They 
are shown in Fig. 7 as a function of the scaled nuclear density. All the data are along the 
smooth curve, whereas the 9 Be slope cannot be explained by such smooth density de- 
pendence. However, if the maximum local density is taken as the abscissa, the situation 
is different as shown in Fig. 8. The curve is plotted to interpolate the shell-model slopes 
except for 9 Be. If the maximum density is calculated by the shell model, the modifica- 
tion slope \dR/dx\ is too large to be expected from the 9 Be density, whereas it is almost 
on the curve if the maximum density is calculated in the AMD with the clustering struc- 
ture. The "anomalous" JLab result on the 9 Be could be explained by the cluster structure 
in the nucleus. We also notice a small difference between the AMD and shell model in 
12 C of Fig. 8, and it is caused by the mixing of a cluster-like configuration in the AMD 
result of 12 C. 

We would like to clarify our viewpoint because it may be confusing to the reader why 
the clustering effects appear in dR/dx of Fig. 8, whereas the effects are small in Fig. 6. 
The nuclear structure functions consist of the mean conventional part and the remaining 
one depending on the maximum local density. The first part is described by the usual 
convolution calculation with the spectral function given by the averaged nuclear density 
distribution, and thus the inhomogeneity of the nuclear density is washed out. The re- 
maining second part is associated with the inhomogeneity of the nuclear density, before 
taking the average of nuclear wave function, given by the nuclear cluster structure. Our 
results imply that the physics mechanism, associated with the high densities created by 
the clusters in 9 Be, could be the origin for the anomalous nuclear-modification slopes 
dR/dx of 9 Be. 



TENSOR STRUCTURE BY QUARK DEGREES OF FREEDOM 



Tensor structure has been investigated in nuclei by hadron degrees of freedom, namely 
by nucleons and mesons. For example, the D-state admixture due to tensor force in 
nucleon-nucleon interactions is the reason for the finite quadrupole moment in the 
deuteron. However, the time has come to investigate the tensor structure by quark 
degrees of freedom. It can be studied by measuring the tensor structure function b\ 
in deep inelastic lepton-deuteron scattering as illustrated in Fig. 9. 



Tensor structure b r for deuteron 



Tensor-structure crisis!? 






only S wave I 1 experiment 

b,= "standard model" b x ^ b "standard model" 

FIGURE 9. Illustration of tensor structure in the deuteron. The structure function b\ 
vanishes if constituents are in the S-wave. A finite b\ can be obtained in the conventional 
deuteron model with the D-state admixture, and it may be considered as the "standard 
model". Nucleon-spin studies imply that such a conventional picture could fail. Then, a 
new field of spin physics will be explored by future experimental measurements. 



It could be predicted in a conventional nuclear-physics approach with the usual con- 
volution description by PDFs convoluted with nucleon-momentum distributions as illus- 
trated in the middle of Fig. 9. However, our experience of nucleon-spin puzzle suggests 
that such a traditional model would not work in spin and orbital angular momenta. There 
is a good possibility that interesting new field of physics will be created by measuring 
the tensor structure functions. 

The structure function b\ is expressed by tensor-polarized distributions bjq as 

bi{x) = + > bit*) = «?(*) - q ' l{x)+ 2 q ' l{x \ (5) 

where i indicates the flavor of a quark, e ; is the charge of the quark, indicates 
an unpolarized-quark distribution in the hadron spin state A, and the distributions are 
defined by the ones per nucleon. For analyzing HERMES data on b\, the distributions 
are expressed as 5jqf v {x) = Srw(x)q?(x) and 5rqf(x) = OCq&rw(x)qf(x). Namely, 
the tensor-polarized distributions are the corresponding unpolarized ones multiplied by 
a weight function 8jw(x) and a constant 0Cq. We consider that a certain fraction of 
the unpolarized distributions are tensor polarized. In expressing the unpolarized PDFs 
in the deuteron, nuclear corrections are neglected, isospin symmetry is used for the 
neutron PDFs in relating them to the proton ones, and flavor-symmetric tensor-polarized 



antiquark distributions are used in drg. The weight function is then parametrized 



&rw(x) = ax b (l —x) c (xq —x), (6) 

where a function with a node is taken by considering a constraint based on a quark- 
parton model: / dxb\(x) = — Jjlirn^o^bW = [21]. Because of this constraint, xq 
is expressed by the other parameters which are determined by a % 2 analysis of the 
HERMES b\ data. We tried two different analyses: 

• Set 1: Tensor-polarized antiquark distributions are terminated (0^=0). 

• Set 2: Finite tensor-polarized antiquark distributions are allowed ((Xq is a parameter). 

The parameters are determined by a x 2 At to the data in Fig. 10, and obtained x 2 values 
are # 2 /d.o.f.=2.83 and 1.57 for set-1 and 2, respectively. As shown in Fig. 10, the set- 
1 analysis is not good enough to explain the HERMES data at small x < 0. 1 without 
the antiquark tensor polarization. The dashed curve is much below the data points, 
whereas the solid curve can describe the data reasonably well. From these analyses, the 
determined tensor-polarized distributions for valence quarks and antiquaries are shown 
in Fig. 1 1 . Of course, the HERMES data are not accurate enough to obtain reliable tensor 
distributions at this stage; however, the obtained distributions can be used for comparing 
them with theoretical calculations and for proposed future experiments such as the ones 
at JLab [22]. 

The b\ sum is expected to vanish in a simple quark-parton model. However, if a 
deviation is found, it suggests a finite tensor-polarized antiquark distributions [21]: 

/ dxb\(x) = ——limtFn(t) + / dx\$8ru(x) +28rd(x) + 8ts(x) + 8ts(x)}, 
J 24f->o ^ 18 7 

/dx If 2 f - 

— [F*(x) -F 2 n (x)] = - / dx[u v (x) - d v [x)} + -J dx[u(x) - d(x)}. (7) 

The b\ sum rule is very similar to the Gottfried one [23], whose violation indicated a 
difference between u and d by the above equation. An interesting point of our analysis 
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FIGURE 10. Two analysis results are shown with 
the HERMES b\ data [14]. The solid (dashed) 
curve indicates the parametrization with (without) 
the antiquark tensor polarization. 



FIGURE 11. Obtained tensor-polarized distribu- 
tions [14]. The solid and dotted (dashed) curves 
indicate the distributions with (without) the anti- 
quark tensor polarization. 



[14] is that a finite sum / dxb\{x) = 0.0058 was obtained to suggest a finite tensor 
polarization for antiquarks. It is interesting to conjecture a possible mechanism of this 
finite antiquark polarization and to test it by an independent experiment. For probing 
antiquark distributions, Drell-Yan process are ideal. In fact, polarized proton-deuteron 
Drell-Yan processes are theoretically formulated in Ref. [24] and they could be measured 
by future hadron facilities such as the J-PARC [25]. The studies of the tensor structure 
function b\ and the tensor-polarized distributions should shed light on a new aspect of 
high-energy hadron spin physics. 
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